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(57) Abstract: A computer (C) is embodied to permit program- 
ming by active objects in multi-task mode, which are represen- 
tative of systems for simulation and comprises a device (D) for 
real world simulation. The device comprises an object simulation 
programme for combined evolution of some at least of the active 
objects, comprising i) objects of state each containing at least one 
item of spatial and/or temporal and/or property data defining a cur- 
rent state, ii) interaction objects each containing the definition of 
at least one of the state objects and at least one function applica- 
ble to at least one of said state objects and defining at each point 
in time the simulated system topology and iii) a simulation man- 
ager which may operate by sequences on a selection of interaction 
objects and to activate each interaction object only once for each 
sequence, according to an order varying in an at least partially ran- 
dom manner from one sequence to the other such as to apply each 
of the functions thereof to the actual state of each state object des- 
ignated thereby for the evolution of the state thereof towards a new 
state. 

(57) Abrege : Un ordinateur (C) est agence de maniere a suppor- 
ter en mode multi-taches une programmation par objets actives, re- 
presentatifs de systemes a simuler, et heberge un dispositif (D) de 
simulation du monde reel. Ce dispositif comprend un logiciel de 
simulation par objets de revolution conjointe de certains au moins 
des objets actives comportant i) des objets d'etat contenant cha- 
cun au moins une donnee d'espace et/ou de temps et/ou au moins 
une donnee de propriete, definissant un etat courant, ii) des objets 
d'interaction contenant chacun la designation d'au moins Tun des 
objets d'etat et d'au moins une fonction applicable a l'un au moins 
de ces objets d'etat, et definissant a chaque instant la topologie du 
systeme simule, et iii) un gestionnaire de simulation capable de 
travailler 
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DISPOSITIF DE SIMULATION DU MONDE REEL PAR TRAITEMENT 
ASYNCHRONE ET CHAOTIQUE 

5 

L' invention concerne la simulation informatique du monde reel, et notamment de T evolution 
temporelle de milieux objets de phenomenes physiques et/ou chimiques et/ou biologiques. 

10 Comme le sait Phomme de Tart, de tres nombreux objets (animes ou inanimes) ou 
substances interagissent avec les milieux (ou media) dans lesquels ils sont places. Or, 
lorsque Phomme con9oit un nouvel objet, comme par exemple un bateau, ou bien lorsqu'il 
tente de comprendre un mecanisme reactionnel, comme par exemple le developpement d'un 
cancer, il lui est generalement impossible de prendre en compte toutes les interactions, ou 

15 du moins les principales, notamment lorsqu' il utilise des maquettes. Par consequent, il arrive 
frequemment que Pobjet final ne procure pas les resultats escomptes, ou que P explication 
avancee ne corresponde pas a la realite, en raison de la complexity et/ou de la combinaison 
des interactions presentes dans le monde reel. Ainsi, un bateau peut s'averer incapable de 
supporter certains etats de mer, alors meme qu'il avait fait Pobjet d'etudes prealables 

2 0 poussees, y compris dans un bassin de carene. 

Afin de mieux anticiper le comportement final d'un objet ou produit, ou de mieux 
comprendre un mecanisme reactionnel, de nombreux industriels et de nombreux chercheurs 
utilisent des dispositifs de simulation informatique du monde reel. 

25 

Ces dispositifs de simulation reposent generalement sur une modelisation des phenomenes 
physiques et/ou chimiques et/ou biologiques sous la forme d' equations differentielles. En 
P absence de solution analytique, la simulation informatique doit s'appuyer sur des methodes 
iteratives, qui passent d'un etat courant du systeme considere a un etat suivant de ce meme 

3 0 systeme. Chaque systeme est generalement decompose en sous-systemes ou mailles dont les 

evolutions respectives sont calculees en parallele, de fa9on synchrone. Par consequent, Petat 
du systeme a P instant T+l resulte de P application en parallele des phenomenes selectionnes 
a Petat de chaque maille a P instant T, sans se preoccuper des consequences que cela peut 
induire entre mailles. En outre, cette technique de simulation ne favorise pas les modifica- 
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tions des parametres de simulation en cours de travail, ni la combinaison de modeles de 
natures differentes, ce qui nuit a sa generality d' application. 

L'invention a done pour but d'ameliorer la situation. 

5 

Elle propose a cet effet un dispositif de simulation du monde reel, propre a etre implante 
dans un ordinateur capable de supporter en mode multi-taches une programmation par objets 
actives, representatifs de systemes a simuler. 

10 Ce dispositif se caracterise par le fait qu'il comprend un logiciel de simulation par objets de 
revolution conjointe de certains au moins des objets actives, comportant : 

- des objets d'etat contenant chacun au moins une donnee d'espace et/ou de temps et/ou au 
moins une donnee de propriete, definissant un etat courant, 

- des objets d'interaction contenant chacun la designation d'au moins Tun des objets d'etat 
15 et d'au moins une fonction applicable a l'un au moins de ces objets d'etat, et definissant a 

chaque instant la topologie du systeme simule, 

- un gestionnaire de simulation capable de travailler par sequences sur une selection d' obj ets 
d'interaction, et d'activer chaque objet d'interaction une unique fois lors de chaque 
sequence, selon un ordre variant de fa<?on au moins partiellement aleatoire d'une sequence 

2 0 a l'autre, de maniere a appliquer chacune de ses fonctions a l'etat courant de chaque objet 
d'etat qu'il designe pour faire evoluer son etat vers un nouvel etat courant. 

En d'autres termes, le dispositif selon l'invention fonctionne selon un mode asynchrone, du 
fait que les etats respectifs des objets actives varient les uns apres les autres au sein de 

2 5 chaque sequence, compte tenu des etats respectifs des autres objets actives, et chaotique, du 

fait que l'ordre de traitement de chaque objet active varie de fagon aleatoire d'une sequence 
a l'autre. II est possible d'activer ou supprimer a tout moment (e'est-a-dire en temps reel) 
un ouplusieurs objets, afin de modifier les conditions de travail et/ou le systeme simule sans 
qu'il faille recommencer integralement la simulation, ce qui confere au dispositif un 

3 0 veritable caractere interactif. 

Le logiciel de simulation peut en outre comprendre des objets d'interaction interne contenant 
la designation d'un seul objet d'etat et d'au moins une fonction applicable a cet objet d'etat, 
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et des objets d' interaction mutuelle contenant la designation d'au moins deux objets d'etat 
et d'au moins une fonction applicable a des donnees de propriete de ces objets d'etat 
designes. 

5 Par ailleurs, certains au moins des obj ets d' etat peuvent comprendre une donnee de propriete 
qui represente une grandeur intensive, et/ou certains au moins des objets d'interaction 
peuvent comporter une fonction faisant intervenir une grandeur extensive ou intensive. 

On peut egalement utiliser des objets d'etat qui pour certains au moins d'entre eux, 
1 0 comprennent des sous-objets d'etat, ainsi qu'eventuellement des sous-objets d'interaction 
operant sur ces sous-objets d'etat. 

En outre, le logiciel de simulation peut comporter des classes d'objets qui definissent des 
structures d'objets d'etat et d'objets d'interaction. Les objets d'etats et d'interaction sont 
1 5 alors derives de ces classes par instanciation. 

De plus, le logiciel de simulation peut comporter un ordonnanceur capable de fonctionner 
soit selon un mode en temps reel* dans lequel il fonctionne selon une frequence choisie, soit 
selon un mode en temps virtuel, dans lequel il fonctionne de fa<?on periodique mais pendant 
2 0 des durees variables d'une periode a l'autre. 

D'autres caracteristiques et avantages de Tinvention apparaitront a Texamen de la description 
detaillee ci-apres, et des dessins annexes, sur lesquels : 

- la figure 1 illustre de fa9on tres schematique, sous la forme de blocs fonctionnels, un 

2 5 ordinateur equipe d'un exemple de realisation d'un dispositif de simulation selon 

Pinvention, et 

- la figure 2 est un organigramme detaillant un exemple de fonctionnement d'un dispositif 
selon P invention. 

3 0 Les dessins annexes pourront non seulement servir a completer Pinvention, mais aussi 

contribuer a sa definition, le cas echeant. 

L'invention concerne un dispositif de simulation du monde reel, interactif. 
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Comme cela est illustre schematiquement et de fa9on fbnctionnelle sur la figure 1, un 
dispositif de simulation selon rinvention D peut-etre installe dans un ordinateur C 
comportant un systeme d' exploitation OS et des moyens de traitement et de calcul CPU 
adaptes aun fonctionnement dans un mode multi-taches, tel que celui offert par Penvironne- 
5 ment oRis decrit notamment dans le document "Systemes multi-agents", pages 499 a 524, 
RSTI - TSI 5 21/2002. Un tel environnement multi-taches est particulierement bien adapte 
a la programmation par objets actives, par exemple en langage C++ ou Java. L' environne- 
ment multi-taches oRis est couple, comme le dispositif D qui est illustre, a un compilateur 
(ici nomme "compilateur de programmation objet"). Par ailleurs, 1 ' environnement oRis peut- 

10 etre couple, comme le dispositif D qui est illustre, a un traducteur en langage C++ (ici 
nomme "interpreteur de programmation objet") afin d'ameliorer son efficacite par 
compilation. Cet interpreteur peut meme etre adapte de maniere a constituer un compilateur 
en ligne dans lequel le code execute est un code compile en ligne et modifiable de fa<?on 
dynamique. Un tel environnement multi-taches, constituant une evolution d' oRis, est connu 

1 5 sous le nom AReVi. 

Le dispositif D comprend un logiciel de simulation de revolution conjointe d' objets actives 
(ici nomme "simulateur general"). Plus precisement, ce logiciel de simulation (ou simulateur 
general) comporte des objets qui peuvent etre decomposes en deux groupes. 

20 

Un premier groupe (nomme "groupe d'objets d'etat") comprend des objets dits "d'etat" 
contenant chacun une ou plusieurs donnees d'espace (X,Y,Z) et/ou de temps (T) et/ou une 
ou plusieurs donnees de propriete (ou parametres, ou encore attributs), qui definissent 
ensemble un etat courant pour l 5 objet d'etat conceme. 

25 

Dans le cas le plus simple, dans lequel il n'existe qu'un seul objet d'etat active, 1'objet 
d'interaction peut ne pas designer 1'objet d'etat, celui-ci etant alors designe de fa<?on 
implicite. 

3 0 Par ailleurs, dans le cas le plus simple, la donnee d'espace-temps d'un objet d'etat peut se 
reduire au seul temps. 
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Certains au moins des objets d' etat comprennent une donnee de propriete qui represente une 
grandeur intensive. 

Dans un premier exemple relatif a un reacteur chimique, un obj et d'etat peut-etre un volume 
5 non localise, et une donnee de propriete peut-etre une concentration d'une substance 
chimique. 

Dans un second exemple relatif a une veine sanguine, un objet d'etat peut-etre une maille 
volumique localisee dans l'espace, une donnee d'espace peut-etre la localisation geome- 
10 trique dans un repere choisi, et une donnee de propriete peut-etre une concentration 
chimique ou la temperature a Pinterieur d'une maille. 

Un second groupe (nomme "groupe d'objets d' interaction") comprend des objets dits 
"d'interaction'' (ou activites) contenant chacun la designation d'au moins l'un des objets 
1 5 d'etat et d'au moins une fonction applicable a l'un au moins de ces objets d'etat designes. 
Le mot "fonction" doit etre compris dans sa definition informatique et non mathematique, 
selon laquelle elle prend quelque chose (n'importe quoi de defini) pour delivrer un resultat. 
Une fonction peut ainsi designer une methode, ou un precede, ou encore un comportement 
a mettre en oeuvre. 

20 

II est important de noter qu'une fonction peut dependre des donnees de propriete d'objets 
d'etat. 

Certains au moins des objets d'interaction comprennent une donnee de propriete qui 

2 5 represente une grandeur intensive ou extensive. 

Dans le premier exemple precite, relatif a un reacteur chimique, un objet d'interaction peut- 
etre une reaction chimique qui manipule les donnees de propriete "concentration(s)" de 
l'objet d'etat "volume", et la reaction chimique peut-etre caracterisee par un parametre 

3 0 "vitesse de reaction" qui peut dependre des donnees de propriete "concentration(s)". 

Dans le second exemple precite, relatif a la veine sanguine, un objet d'interaction peut-etre 
un "diffiiseur thermique" qui equilibre les donnees de propriete de temperature entre deux 
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mailles adjacentes, et le diffuseur thermique peut-etre caracterise par le parametre 
"conductivity thermique". 

Chaque objet d'interaction (ou activite) est done associe a au moins un objet d'etat par une 
fonction, et un objet d'etat peut etre associe a plusieurs objets d'interaction (ou activites). 
Une fonction ne modifie que les parametres du ou des objets d'etat auxquels elle est 
associee, sans modifier les parametres des autres objets d'etat et objets d' interaction. Par 
ailleurs, les caracteristiques des objets d'etat peuvent permettre de choisir la fonction 
associee aux objets d'interaction ou d'en modifier les caracteristiques. En outre, une 
fonction peut-etre modifiee de fa(?on temporaire ou permanente en fonction des parametres 
d'un objet d'etat associe. 

Comme on le verra plus loin, un objet d'interaction peut constituer une interface definie, 
generalement avec une orientation, entre deux milieux dits "source" et "cible" (eventuelle- 
15 ment confondus) constituant des objets d'etat. L'interface fait alors le lien entre les deux 
milieux, qui sont en liaison a un instant donne (notion de synchronicite forte), et modifie 
leurs parametres respectifs (par exemple, l'interface peut etre une difference de densite entre 
deux milieux qui provoque une diffusion moleculaire a une certaine vitesse). 

Tout type d' interface peut etre envisage, et notamment des interfaces de transport (par 
exemple de molecules, de chaleur, de courant, et analogues), et des interfaces de relaxation 
(par exemple de reactions chimiques ou nucleaires). Par consequent, un milieu peut ne pas 
avoir d' extension spatiale, celle-ci etant induite, tout comme les eventuelles vitesses, par une 
ou plusieurs interfaces. Par ailleurs, une activite appliquee a un milieu peut modifier le 
milieu, comme par exemple des phenomenes de relaxation, ce qui correspond a une auto- 
activation interne au milieu. 

Le logiciel de simulation (ou simulateur general) comprend egalement un gestionnaire de 
simulation couple aux groupes d'objets d'etat et d'objets d'interaction et agence de maniere 
3 0 a creer son propre sequenceur ou ordonnanceur (ou en anglais "scheduler") afin de travailler 
de fa9on sequentielle sur une selection d'objets d'interaction dudit groupe d'objets 
d'interaction. Plus precisement, le gestionnaire de simulation est charge d'activer une unique 
fois lors de chaque sequence, sous le controle du sequenceur (ou scheduler) qu'il cree pour 
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P occasion, chaque objet d' interaction selectionne, selon un ordre qui varie de fa9on au 
moins partiellement aleatoire d'une sequence a P autre, afin d'appliquer chacune de ses 
fonctions a Petat courant de chaque objet d'etat qu'il designe de maniere a faire evoluer son 
etat vers un nouvel etat courant. 

5 

En d' autres termes, comme cela est iliustre a titre d' exemple sur 1 ' organi gramme de la figure 
2, Putilisateur choisit tout d'abord, dans une etape 10 5 parmi les premier et second groupes 
d'objets, un ou plusieurs objets d'etat et un ou plusieurs objets d'interaction qui se 
rapportent a ce ou ces objets d'etat choisis, afin que le dispositif D simule P evolution spatio- 
1 0 temporelle du systeme represents par lesdits objets choisis. Cela "pre-active" chaque objet 
d' interaction choisi au sein du logiciel de simulation. 

Puis, dans une etape 20, le gestionnaire de simulation initialise un compteur de sequences 
en pla9ant la valeur n du compteur a 1 , et cree une liste d'objets d'etat et d' objets d'interac- 
15 tion. La premiere sequence (n=l) commence par une etape 30 dans laquelle le gestionnaire 
de simulation selectionne de fa9on aleatoire, au sein de la liste d'objets d' interaction, Pun 
des objets d' interaction choisis, ce qui Pactive momentanement. Chaque fonction de Pobjet 
d'interaction selectionne et active est alors appliquee a Petat courant de chaque objet d'etat 
de la liste d'objets d'etat avec lequel il est en relation, ce qui modifie eventuellement son 

2 0 etat en cours. L'activation et Papplication (ou execution) de chaque fonction a chaque objet 

d'etat choisi constituent Petape 40. 

Puis, dans une etape 50, le gestionnaire de simulation supprime de la liste d'objets d'interac- 
tion de la sequence en cours Pobjet d'interaction qu'il vient d'appliquer. 

25 

Ensuite, dans une etape 60, le gestionnaire de simulation effectue un test destine a 
determiner s'il reste d'autres objets d'interaction a appliquer dans la liste d'objets d'interac- 
tion de la sequence en cours. 

3 0 Si la liste d'objets d'interaction de la sequence en cours n'est pas vide, le gestionnaire de 

simulation retourne a Petape 30 afin de proceder a la selection de fa9on aleatoire d'un objet 
d'interaction restant. Comme indique ci-avant, il active alors ce nouvel objet d'interaction 
selectionne et applique chacune de ses fonctions a Petat courant de chaque objet d'etat, 
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eventuellement modifie par 1'activation de Pobjet d'interaction precedent (lequel ne peut 
plus etre utilise dans la sequence en cours). Le gestionnaire de simulation reproduit ces 
operations (selection-activation, application et mise(s) a jour) autant de fois qu'il y a 
d'objets d'interaction selectionnes dans sa liste d'objets d' interaction, de sorte que chaque 
fonction de chaque objet d' interaction soit appliquee une unique fois k chaque objet d'etat 
active auquel il est associe. Une fois celatermine, la premiere sequence (n=l), correspondant 
a Tinstant TK), est achevee. 

Lorsque la liste d'objets d'interaction de la sequence en cours est vide, le gestionnaire de 
simulation incremente d'une unite, dans une etape 70, la valeur en cours n du compteur de 
sequences. Bien entendu, cette etape 70 comporte un test sur le nombre de sequences a 
effectuer. Si le nombre de sequences effectuees est egal au nombre maximal prevu, le 
gestionnaire de simulation met fin a la simulation. En revanche, s'il reste au moins une 
sequence a effectuer, le gestionnaire de simulation retourne a l'etape 20 pour effectuer une 
nouvelle sequence correspondant a un instant T+l, T+2, T+n. II reitere alors a chaque 
nouvelle sequence les operations precitees. 

La duree de la simulation, et done le nombre maximal de sequences effectuees par le 
gestionnaire de simulation, depend de V application concernee, ou duparametrage choisi par 
l'utilisateur compte tenu de Papplication. Mais, la simulation peut etre interrompue a chaque 
instant par V utilisateur a Taide d'une instruction d 5 arret transmise au logiciel de simulation 
grace a une interface homme/machine de Tordinateur C. II est important de noter qu'une 
simulation interrompue a la demande d'un utilisateur peut-etre reprise ulterieurement. 

En raison de ce mode de fonctionnement sequentiel et chaotique du logiciel de simulation, 
Tutilisateur peut a chaque instant intervenir dans une simulation, soit sous la forme d'un 
"avatar" pour interagir Iui-meme avec 1'objet de la simulation, par exemple le milieu, soit 
pour ajouter a, ou supprimer de, sa selection un ou plusieurs objets d'etat et/ou un ou 
plusieurs objets d'interaction. L'utilisateur peut egalement decider de modifier au moins 
partiellement la definition (ou structure) d'un ou plusieurs objets d'etat ou d'interaction. 
Cela confere au logiciel de simulation une grande interactivite. 
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Trois exemples permettant de mieux comprendre ce que representent les objets d'etat et 
d' interactions sont donnes ci-apres. Pour eviter d'alourdir inutilement la description, les 
exemples sont simplifies et reduits a ce qui est necessaire pour permettre a l'homme du 
metier de comprendre comment on peut mettre en oeuvre 1' invention. 

5 

Un premier exemple conceme la simulation de la cinetique chimique au sein d'un reacteur 
chimique. Le reacteur chimique constitue un objet d'etat representant un milieu comportant 
des substances chimiques, et au sein duquel peuvent survenir N reactions chimiques 
constituant autant d' objets d' interaction. 

10 

L' objet d'etat reacteur chimique est par exemple associe a des donnees de propriete 
representant les concentrations (C1(T), C2(T), Cn(T)) a l'instant T des n substances 
chimiques presentes initialement dans le reacteur. Ces concentrations definissent ici l'etat 
(a Pinstant T) de l'objet d'etat reacteur chimique, dont on souhaite simuler revolution 
15 temporelle. 

Chaque objet d'interaction reaction chimique est par exemple caracterise par des donnees 
de propriete representant, d' une premiere part, la vitesse de la reaction V = f(C 1 , C2, . . . , Cn), 
d'une deuxieme part, les reactifs concernes par la reaction (Rl , R2, Rn), d'une troisieme 
2 0 part, les produits resultants de la reaction (PI, P2, Pm), et d'une quatrieme part les 
coefficients stoechiometriques definissant les proportions initiales des reactifs a T==0. En 
outre, chaque objet d'interaction reaction chimique est par exemple associe a une fonction 
representant le comportement de la reaction : dans un premier temps on lit les concentrations 
(C1(T), C2(T), Cn(T)) en cours des reactifs presents a l'instant T dans le reacteur 

2 5 chimique, dans un deuxieme temps on calcule la vitesse de reaction V qui depend des 

concentrations des reactifs presents dans le reacteur chimique, et dans un troisieme temps 
on modifie les concentrations des reactifs et des produits apres un temps dT a la vitesse de 
reaction V. 

3 0 Ici chaque objet d'interaction n'agit que sur un unique objet d'etat, de sorte qu'il constitue 

un objet d'interaction interne. 



WO 2005/081141 PCT/FR2005/000124 

10 

Dans ce premier exemple, l'utilisateur va done choisir son objet d'etat (e'est-a-dire la 
composition initiale de son reacteur chimique) et sesN objets d'interaction (e'est-a-dire les 
N reactions chimiques qui sont impliquees du fait de la composition choisie). Le 
gestionnaire de simulation peut alors debuter son traitement sequentiel. 

5 

Comme indique precedemment, il debute une premiere sequence en selectionnant de fa<?on 
aleatoire Pun des N objets d'interaction choisis. Get objet d'interaction est active afin qu'il 
applique sa fonction (comportement de la reaction) a l'etat courant de l'unique objet d'etat, 
defini ici par les concentrations (C1(T), C2(T), Cn(T)) en cours des reactifs (lors de la 

10 premiere sequence T=0). Le gestionnaire met a jour les concentrations (C1(T), C2(T), 

Cn(T)) des reactifs. Puis, il selectionne de fa9on aleatoire l'un des N-l objets d'interaction 
restants, afin de l'activer et d'appliquer sa fonction au nouvel etat courant de l'unique objet 
d'etat. II reproduit ces operations N fois de sorte que chaque fonction de chaque objet 
d'interaction soit appliquee une unique fois a l'objet d'etat. La premiere sequence est alors 

1 5 terminee. 

Le gestionnaire debute alors une nouvelle sequence, correspondant a Finstant T+l, si cela 
s'avere neeessaire, en reiterant les operations precitees. La simulation prend fin lorsqu'il ne 
reste plus de reactifs pour produire des produits, ou bien lorsque l'utilisateur adresse une 
2 0 instruction d' arret au logiciel de simulation. 

Un deuxieme exemple concerne la simulation de la diffusion moleculaire de M diffuseurs 
au sein d'un espace decoupe en K mailles. 

2 5 Chacune des K mailles constitue un objet d'etat caracterise, par exemple, par une donnee 
de position (ou localisation topographique ou geographique dans un repere 3D) et par des 
donnees de propriete representant les concentrations (C1(T), C2(T), Cn(T)) a Pinstant 
T de n substances chimiques (SI, S2, Sn) presentes dans la maille. Ces concentrations 
definissent ici l'etat (a l'instant T) de la maille. 

30 

Chacun des M diffuseurs constitue un objet d'interaction caracterise, par exemple, par des 
donnees de propriete (ou parametres) representant, d'une premiere part, la substance diffusee 
Si, d'une deuxieme part, la vitesse de diffusion tridimensionnelle V = (Vx, Vy, Vz) de la 
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substance Si, et d'une troisieme part, les deux mailles entre lesquelles s'effectue la diffusion. 
En outre, chaque diffuseur (objet d' interaction) est par exemple associe a une fonction 
representant son comportement : dans un premier temps on lit la concentration Ci(T) de la 
substance Si a diffuser dans chacune des deux mailles, dans un deuxieme temps on calcule 
5 les quantites a diffuser (par exemple en utilisant la loi de Fick generalisee), et dans un 
troisieme temps on modifie la concentration du diffuseur Si dans chacune des deux mailles 
apres un temps dT a la vitesse de diffusion V. Les donnees de propriete representent ici les 
attributs (ou parametres) de la fonction (comportement du diffuseur). 

10 Ici chaque objet d' interaction agit sur deux objets d'etat, de sorte qu'il constitue un objet 
d' interaction mutuelle. Par ailleurs, chaque objet d' interaction presente ici une structure (ou 
definition) classique reposant sur au moins une fonction designant au moins un objet d'etat, 
ainsi qu' eventuellement des parametres ou regle(s), ou encore loi(s), lies a la fonction et 
constituant des donnees de propriete. Mais, leur structure peut etre plus complexe, par 

1 5 exemple lorsqu'elle comporte un ou plusieurs sous-objets d'interaction. En outre, les mailles 
sont ici de type tridimensionnel (3D), mais dans certaines applications elles peuvent etre de 
type bidimensionnel (2D) ? voire meme unidimensionnel (ID). 

Dans ce deuxieme exemple, Futilisateur va done choisir ses K objets d'etat (e'est-a-dire les 
2 0 K mailles de diffusion) et ses M objets d'interaction (e'est-a-dire les M diffuseurs). Le 
gestionnaire de simulation peut alors debuter son traitement sequentiel. 

II debute une premiere sequence en selectionnant de fafon aleatoire Tun des M objets 
d'interaction choisis. Cet objet d'interaction est active afin qu'il applique sa fonction 

2 5 (comportement du diffuseur), compte tenu de ses propres attributs, a Pet at courant des K 

objets d'etat (mailles), defmi ici par les concentrations (C1(T), C2(T), Cn(T)) en cours 
des substances chimiques (SI, S2, Sn) dans les differentes mailles (lors de la premiere 
sequence T=0). Le gestionnaire met a jour les concentrations (C1(T), C2(T), Cn(T)) des 
substances (SI, S2, Sn) dans chacune des K mailles. Puis, il selectionne de fa9on 

3 0 aleatoire l'un des N-l objets d'interaction restants, afin de l'activer et d'appliquer sa 

fonction au nouvel etat courant de chaque objet d'etat. II reproduit ces operations N fois de 
sorte que chaque fonction de chaque objet d'interaction soit appliquee une unique fois a 
chaque objet d'etat. La premiere sequence est alors terminee. 
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Le gestionnaire debute alors une nouvelle sequence, correspondant a l'instant T+l, si cela 
s'avere necessaire, en reiterant les operations precitees. La simulation prend fin lorsque les 
concentrations (C1(T), C2(T), Cn(T)) des substances chimiques (SI, S2, Sn) sont 
respectivement identiques dans chacune des K mailles, ou bien lorsque Putilisateur adresse 
5 une instruction d' arret au logiciel de simulation. 



Dans un troisieme exemple., on peut combiner les deux exemples precedents de maniere a 
simuler revolution de la cinetique chimique et de la diffusion moleculaire au sein d'une 
meme veine. 

10 

Par exemple, on choisit K objets d'etat constituant N reacteurs chimiques localises dans 
1 ' espace 3D. Chaque reacteur chimique est par exemple associe a une donnee de position (ou 
localisation geographique dans un repere 3D) et a des donnees de propriete representant les 
concentrations (C1(T), C2(T), Cn(T)) a Finstant T de n substances chimiques (SI, S2, 
15 Sn) presentes dans la maille. Ces concentrations definissent ici Petat (a Tinstant T) de la 

maille. 

Par ailleurs, on choisit par exemple N+M objets d'interaction constituant N reactions 
chimiques et M dififuseurs. 

20 

Chaque reaction chimique est par exemple associee a des donnees de propriete representant, 
d'une premiere part, la vitesse de la reaction V = f(Cl , C2, Gn), d'une deuxieme part, les 
reactifs concemes par la reaction (Rl, R2, Rn), a T==0, d'une troisieme part, les produits 
resultants de la reaction (PI, P2, Pm), et d'une quatrieme part les coefficients 

2 5 stoechiometriques definissant les proportions initiales des reactifs a T=0. Chaque reaction 

chimique est egalement, par exemple, associe a une fonction representant le comportement 
de la reaction : dans un premier temps on lit les concentrations (C1(T), C2(T), Cn(T)) en 
cours des reactifs presents a 1' instant T dans le reacteur chimique, dans un deuxieme temps 
on calcule la vitesse de reaction V qui depend des concentrations des reactifs presents dans 

3 0 le reacteur chimique, et dans un troisieme temps on modifie les concentrations des reactifs 

et des produits apres un temps dT a la vitesse de reaction V. 
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Chaque diffuseur est par exemple associe a des donnees de propriete representant, d'une 
premiere part, la substance diffusee Si, d'une deuxieme part, la vitesse de diffusion 
tridimensionnelle V = (Vx, Vy, Vz) de la substance Si, et d'une troisieme part, les deux 
mailles entre lesquelles s'effectue la diffusion. Chaque diffuseur est egalement, par exemple, 
5 associe a une fonction representant son comportement : dans un premier temps on lit la 
concentration Ci(T) de la substance Si a diffuser dans chacune des deux mailles, dans un 
deuxieme temps on calcule les quantites a diffuser (par exemple en utilisant la loi de Fick 
generalisee), et dans un troisieme temps on modifie la concentration du diffuseur Si dans 
chacune des deux mailles apres un temps dT a la vitesse de diffusion V. Les donnees de 
10 propriete represented ici les attributs (ou parametres) de la fonction (comportement du 
diffuseur). 

Le fonctionnement du dispositif est alors identique a celui presente dans les examples 
precedents, mais cette fois les concentrations des differentes substances evoluent non 
15 seulement en fonction des reactions chimiques, mais egalement en fonction des mailles ou 
elles s'effectuent. Par consequent, chaque sequence comporte N+M tirages aleatoires 
permettant d'appliquer successivement les fonctions des N+M objets d'interaction aux K 
reacteurs chimiques. 

2 0 Bien entendu, 1'invention peut s'appliquer a des applications beaucoup plus complexes que 
celles presentees ci-dessus a titre d' exemple illustratif. Elle s'applique egalement a des 
applications plus simples dans lesquelles le logiciel n'intervient que sur un unique objet 
d'etat active, a l'aide d'un unique objet d'interaction interne. Par ailleurs, dans certaines 
applications Tun au moins des objets d'etat peut presenter une structure (ou definition) 

2 5 complexe reposant sur un ou plusieurs sous- objets d'etat (presentant la structure classique 

d'un objet d'etat), eventuellement associe a un ou plusieurs sous-objets d'interaction 
(presentant la structure classique d'un objet d'interaction). 

II est maintenant fait reference au tableau de Tannexe XI, pour une description, au niveau 

3 0 du code, d'un exemple de mise en oeuvre de Tinvention faisant intervenir un objet d'etat 

appele "milieu", et plus precisement "classe milieu", et un objet d'interaction appele 
"interface", et plus precisement "classe interface". Le but de cette application est d'associer 
une indication de couleur a une indication de densite au sein du milieu. 
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Le code n'apparait que dans la seconde colonne du tableau, dont la premiere colonne ne 
comporte que des reperes alphanumeriques ne faisant pas partie du code. 

L'annexe XI commence, en A, par des declarations "#include" qu'il est inutile de detailler 
5 du fait qu'elles sont bien connues de Phomme du metier. Ensuite, en B, est declaree une 
classe interface, associee a la classe milieu traitee en detail en C. Apres Pen-tete de classe 
CI, la classe milieu comprend la declaration de methodes specifiques C20 a C30, puis en 
C4 une fonction (ou methode) "activate", et enfin des variables protegees enoncees en C5. 
Dans l 5 annexe XI, les declarations de methodes (ou fonctions) sont allegees des definitions 
10 de type/parametres, puisque Phomme de metier retrouvera celles-ci dans Penonce detaille 
de chaque methode. 

Parmi les variables protegees indiquees en C5, on peut noter : 

- des "_interfaces", qui sont regroupees au sein d'un tableau d'un type particulier dit 
1 5 StlVector (classe standard C++), et 

- une variable "Jf \ qui est de type ArRef (classe propre a Penvironnement AReVi). 

Vient ensuite le detail de la classe interface (par souci d'homogeneite elle reste designee par 
B). Elle comprend une methode de definition et une methode de construction, suivies de 
2 0 declaration de methodes specifiques B20 a B3 1 et de la methode activate B4. Enfin sont 
definies des variables protegees B5 ou Ton note a nouveau en particulier : 

- une "jsource" qui est de type milieu, et 

- une " Jarget" qui est egalement de type milieu. 

2 5 Les definitions des methodes de la classe milieu interviennent ensuite (en C 1 1 -C 1 4 et C20- 

C3 0) en grand detail. Les rubriques C 1 1 a C 1 3 definissent le constructeur de la classe milieu. 
En CI 1 sont initialisees les variables protegees de la classe milieu. 

En C12 apparait une definition d'une forme 3D dite "Vrml 11 . Le code Vrml permet de faire 

3 0 de la visualisation 3D en ligne. Cette forme est definie par une chaine de caracteres 

conformement a la syntaxe Vrml decrite par exemple a P adresse http://www.vrml.org. Cette 
forme est ensuite visualisee et mise en couleurs. 
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La rubrique C 1 3 associe au milieu une activite destinee a declencher la methode activate de 
larubrique CI 4. 

La rubrique C14 est un destructeur de Pobjet d'etat milieu. 

5 

Les rubriques C20 a C30 incluse donnent le detail de la construction des methodes 
specifiques deja evoquees sous la meme identite pour la classe milieu. 

La methode activate a appliquer au milieu est definie en C4. Elle permet de relier une 
1 0 indication de couleur a Tindication de densite, et plus precisement de changer la couleur par 
tirage aleatoire. En d'autres termes ? on modifie les attributs rouge, vert, bleu du milieu 
chaque fois qu 5 activate est appliquee (ou activee). 

Interviennent ensuite des definitions pour la classe interface, avec en B10-B1 1 la definition 
15 du constructeur. Plus precisement, en BIO on initialise les variables protegees de la classe 
interface, tandis qu'en Bll on associe a P interface une activite destinee a appliquer la 
methode activate definie en B4. 

La rubrique B12 est un destructeur de Pobjet d'interaction interface. 

20 

Viennent ensuite les details des methodes specifiques B20 a B31 deja citees. 

La methode activate de ¥ interface est definie en B4. II s'agit de calculer une diffusivite en 
fonction d'une densite de milieu source et d ! une densite de milieu cible, afin de propager 

2 5 Tecart entre les deux densites, si du moins cet ecart est significatif. L'ecart note "delta" est 

calcule comme Implication d'une fonction diffusivite a la densite de milieu source rhos, 
diminuee de la densite de milieu cible rhot. Pour ce faire, on demande d'abord au milieu 
source et au milieu cible quelles sont leurs densites respectives. Puis, on calcule la valeur 
absolue des differences de couleurs. Si cette valeur absolue est superieure a une valeur 

3 0 choisie (ici egale a 0,01), on augmente la densite la plus faible et on diminue la densite la 

plus forte. Puis, on effectue la mise a jour des valeurs a une vitesse choisie. 
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Les rubriques B60-B63 permettent de definir le type des fenetres de visualisation 3D qui 
seront utilisees par P application. 

C'est alors qu'intervient un element important pour Tinvention qui est Pordonnanceur ou 
5 sequenceur (ou en anglais "scheduler"). Les rubriques S10, S20, S21, S22, S30, S31 et S41 
permettent d' initialiser la simulation. 

En l'espece, la rubrique S10 initialise un scheduler en temps virtuel, c'est-a-dire que son 
fonctionnement n'est pas astreint a respecter le temps reel, mais chacune de ses iterations 
1 0 represente logiquement, et non physiquement, une duree de une milliseconde (1 ms). Bien 
entendu, le scheduler peut egalement fonctionner en temps reel. Dans ce cas, chacune de ses 
iterations dure physiquement une periode choisie. 

La rubrique S20 permet d'initialiser les dimensions du milieu grace aux arguments de la 
1 5 ligne de commande. La rubrique S2 1 permet de creer et d' initialiser les milieux. La rubrique 
S22 permet de creer et d'initialiser les interfaces. 

Les rubriques S3 1 et S32 realisent un affichage auxiliaire, qui se fait ici par la voie speciale 
de communication avec l'ecran dit "flot d'erreur". 

20 

La rubrique S41 affiche la scene a traiter, dans des conditions de perspective choisies. 

Enfin, la rubrique S50 est le programme principal main qui initialise le fonctionnement du 
systeme. Plus precisement, ce programme cree tout d'abord un systeme 3D, puis il Pappelle 
25 pour qu'il active la methode donnee a la rubrique S10 afm qu'elle cree le scheduler (ou 
sequenceur), les milieux et Pinterface. 

On comprendra que la structure de classes, presentee ci-avant, peut etre resumee de la 
maniere indiquee schematiquement dans P annexe X2. 

30 

La transition entre la representation de P annexe X2 et le code detaille de Pannexe XI est 
consideree comme accessible a Phomme du metier, des lors qu'un exemple en a ete donne. 
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L'invention trouve de tres nombreuses applications dans de nombreux domaines techniques, 
et notamment dans les domaines de la chimie, de la pharmacie, de la physique, de 
Paeronautique, de P architecture, en particulier navale (par exemple pour 1' etude comporte- 
mental d'un bateau ou d'une plate-forme offshore, en remplacement et/ou en complement 
5 des bassins de carene), de la medecine, notamment dans le cadre de 1' etude du developpe- 
ment et du traitement de certaines maladies (par exemple les cancers) ou de mecanismes 
reactionnels (par exemple l'activation de l'insuline), de Pergonomie, notamment pour la 
realisation d'equipements specifiquement adaptes a des personnes handicapees, et dans le 
domaine de la circulation routiere. 

.0 

L'invention ne se limite pas aux modes de realisation de dispositif de simulation decrits ci- 
avant, seulement a titre d'exemple, mais elle englobe toutes les variantes que pourra 
envisager l'homme de Tart dans le cadre des revendications ci-apres. 
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Annexe XI 



A 


#include statements 




B 


class Interface; 








class Milieu : public Object3D 

{ 

public: 
AR_CLASS(Milieu) 
AR_CONSTRUCTOR(Milieu) 












C20 
C21 


setDensity() 
getDensityQ 




C22 
C23 


setColorQ 
getColorO 




C24 
C25 


addSourceQ; 
remove Source(); 




C26 
C27 
C28 


addlnterface() 

removelnterfaceO 

getNbInterfaces() 




C29 
C30 


getlnterfaceO 
accesslnterface() 




C4 


activate(ArRef<Activity> act, double dt); 




C5 


protected: 

StlVector<ArRef<Interface> > interfaces; 
ArRef<VrmlMaterialInteractor> _it; 
double _red,_green,_blue; 
double ^density; 
bool _changeDensity; 
double isSource; 

>; 




Bl 


class Interface : public ArObject 

{ 

public: 
AR_CLASS(Interface) 
AR_CONSTRUCTOR(Interface) 
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B20 


setVelocity() 




B21 


eetVelocitvO 






o /a^Tl i "fin i o 1 1 / i +i / ( \ 

bcuL/iiiusivi ly \) 




B23 


getDiflusivityQ 




"DO/I 


setSource() 




dZj 


getoource(j 




B26 


access Source() 




JtSZ / 


set i arget(j 




B28 


eetTareetO 




B29 


accessTarget() 




B30 


ff etOth erFn d H 




B31 


accessOtherKn df ^ 




B4 


activate(ArRef<Activity> act, double dt); 




B5 


protected: 






ArRef<Milieu> _sourcej 






ArRef<Milieu> target* 






double velocity diffusivitv* 






AK__CLiAaJS_DEF(Miueu,0 bj ect3D) 




Cll 


Milieu::Milieu(ArCW & arCW) 






: Object3D(arCW), 






_interfaces(), 






Jt(new_VrmlMaterialInteractor()), 






_red(1.0), 






_green(1.0), 






__blue(1.0), 






_density(0), 






changeDensity(true) ? 






_isSource(false) 
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C12 


{ 

ArRef<VnnlShape3D> sh=new_VrmIShape3D0; 
sh->parseString( M \ 
Shape { \ 
appearance Appearance { . \ 
material DEF MAT Material { \ 
diffiiseColor L0 1.0 1.0 \ 
transparency 0.5 \ 
} \ ' 
} \ 
geometry Box { \ 
size 0.9 0.9 0.9 \ 

! * ^ 
}"); 




C13 


ArRef<Activiry> act=new__Activiry(0.1); 
act->setBehavior(thisRefO,&Milieu::activate); 

> 




C14 


Milieu: :~Milieu(void) 

{ 

} 




C20 


void 

Milieu: :setDensity(double density) 
{ 

_density = density; 
_changeDensity = true; 
} 






void 

Milieu: :getDensity(double & density) 
{ 

density = density; 
} 
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C22 


void 

Milieu::setColor(double red, 
double green, 
double blue) 

{ 

_red=red; 

_greerp=green; 

__blue=blue; 

_it->writeDiffliseColorField(_red,_green 5 _blue); 
} 




C23 


void 

Milieu: :getColor(double & redOut, 
double & greenOut, 
double & blueOut) const 

{ 

redOut= red; 
greenOut= __green; 
blueOut=_blue; 
} 




C24 


void 

Milieu: :addSource(void) 
{ 

isSource - true; 

} 




C25 


void 

lVllilC-Ul . .1 CllltJ VCOUUi V'C^VUIU J 

{ 

_isSource = false; 

> 






void 




C26 


Milieu: :addInterface(ArRef<Interface> interface) 
{ 

interfaces.push_back(interface); 

} 
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C27 


void 

Milieu: :removeInterface(ArRef<Interface> interface) 
{ 

for(unsigned int i= interfaces.size();i--;) 
{ 

if(_interfaces[i]— interface) 
{ 

StlVectorFastEraseC_interfaces 5 i); 
break; 

> 

} 

> 




C28 


unsigned hit 

Milieu::getNbInterfaces(void) const 
{ 

return(_interfaces.size()); 

> 




C29 


ArConstRef<Interface> 

Milieu: :getlnterface(unsigned int index) const 

{ 

return(_interfaces [index]); 

> 




C30 


ArRef<Interface> 

Milieu: :accesslnterface(unsigned int index) 

{ 

return(_interfaces[index]); 
} 
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C4 


bool 

Milieu: :activate(ArRef<Activity> act, 
double /* dt */) 

{ 

if(_isSource) { density — 1.0; } 
if(_changeDensity) 
{ 

setColor( 1.0,1 .0-_density 5 1 .0-_density); 

/* 

ifCdensity < 0.1) { setColor(l 5 l 5 l); } 
else { setColor(l 5 0 5 0); } 

*/ 

changeDensity = false; 

} 

return(true); 
} 












BIO 


AR_CLASS_DEF(Interface ? ArObject) 

Interface: :Interface(ArCW & arCW) 
: ArObject(arCW), 

source(), 
J:arget() 9 

velocity(O), 
_dif¥usivity(0) 




Bll 


{ 

ArRef<Activity> act=new_Activity(0. 1); 
act->setBehavior(thisRefO,&Interfece::activate); 

} 




B12 


Interface::~Interface(void) 

{ 

} 




B20 


void 

Interface: :setVelocity(double velocity) 
{ 

velocity = velocity; 

} 
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B21 


void 

Interface::getVelocity(double & velocity) 

{ 

velocity = _velocity; 

> 




B22 


void 

Interface: :setDifrusivity(double diffusivity) 
{ 

_difrusivity = diffusivity; 
} 




B23 


void 

Interface: :getDif¥usivity(double & diffusivity) 
{ 

diffusivity = diffusivity; 
} 




B24 


void 

Interface: :setSource(ArRef<Milieu> source) 

{ 

if(_source) 

{ 

_source->removeInterface(thisRef()); 
} 

source=source; 
if(_source) 
{ 

source->addInterface(thisRef()); 

} 

> 




B25 


ArConstRef<Milieu> 
Interface: :getSource(void) const 
{ 

return(_source); 

} 




B26 


ArRef<Milieu> 

Interface: :accessSource(void) 

{ 

return(_source); 
} 
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B27 


void 

Interface: :setTarget(ArRef<Milieu> target) 
{ 

if(_target) 
{ 

_target->removeInterface(thisRef()) ; 
} 

_target=target; 
if(_target) 

{ 

_target->addInterface(thisRefO); 
> 

} 




B28 


ArConstRef<Milieu> 
Interface: :getTarget(void) const 
{ 

return(_target); 
} 




B29 


ArRef<Milieu> 

Interface: :accessTarget(void) 

{ 

returnC_target); 
} 




B30 


ArConstRef<Milieu> 

Interface: :getOtherEnd(ArConstRef<Milieu> m) const 
{ 

reton(m=_source ? _target : _source); 

> 




B31 


ArRef<Milieu> 

Interface: :accessOtherEnd(ArConstRef<Milieu> m) 
{ 

retura(m==_source ? _target : source); 
} 
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B4 


bool 

Interface: :activate(ArRef<Activity> /* act */, 
double /*dt*/) 

{ 

double rhos 3 rhot; 
_source->getDensity(rhos); 

target->getDensity(rhot); 
bool change=false; 
if(fabs(rhot-rhos)>0.0 1) 

{ 

double delta=_diffusivity*(rhos-rhot); 

rhos -= delta; 

rhot += delta; 

change=true; 

} 

if(change) 
{ 

_source->setDensity(rhos) ; 

_target->setDensity(rhot); 

} 

return(true); 
} 




B60 


class MyViewer : public ExaminerSD 
{ 

public: 

AR^CLASS(MyViewer) 

AR_CONSTRUCTOR(MyViewer) 
protected: 

virtual void onKeyPress(const Viewer3D::KeyPressEvent & evt); 

}; 




B61 


AR_CLASS_DEF(MyViewer ? Examiner3D) 

MyViewer: :MyViewer(ArCW & arCW) 

: Examiner3D(arCW) 

{ 

} 





WO 2005/081141 



PCT/FR2005/000124 



27 



B62 


My Viewer : :~My Viewer(void) 

{ 
} 




B63 


void 

My Viewer: :_onKeyPress(const Viewer3D: :KeyPressEvent & evt) 
{ 

Examiner3D : :_pnKeyPress(evt) ; 
if((evt.key="+")||(evt.key="-")) 
{ 

double dist,yaw,pitch; 

computeCursorDirection(evt.xMouse,evt.yMouse,yaw J pitch); 

ArRef<Base3D> ray=new_Base3D(); 

ray->setLocation(thisRef()); 

ray->yaw(yaw); 

ray->pitch(pitch); 

ArRef<Object3D> found; 

getScene()->firstRayMersection(ray ? 3 3 true,Milieu::CLASS() ? 
found 9 dist,thisRefQ); 

if(found) 

{ 

ArRef<Milieu>m(static_cast<Miiieu*>(found.c_ptr())); 

cerr « "Found '* « m->getAlias() « " M « evt.key « endl; 

if(evt.key==="+") { m->addSource(); } 

if(evt.key="-") { m->removeSourceO; } 

} 

} 

} 




S10 


ArRef<Scheduler> 
simulationInit(void) 

{ 

ArRef<Scheduler> sched^ew_VirtualTimeScheduler(l e-3); 
unsigned int nbX=40,nbY=20 s nbZ=2; 
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S20 


if(ArSystem: :getCommandLine().size()> 1 ) 
{ 

strToUint(nbX,ArSystem::getCommandLine()[l]); 
} 

if(ArSystem::getCommandLine().size()>2) 
{ 

strToUint(nbY,ArSystem::getCommandLine()[2]); 
} 

if(ArSystem: : getCommandLine() . s ize()>3 ) 
{ 

strToUint(nbZ 5 ArSystem::getCommandLine()[3]); 
} 




S21 


unsigned int x,y,z; 

StlVector<StlVector<StlVector<ArRef<Milieu> > > > t; 
for(x=0;x<nbX;x++) 
{ 

StlVector<StlVector<ArRef<Milieu> > > ex; 
t.push_back(ex); 
for(y=0;y<nbY;y++) 
{ 

StlVectoKArRef<Milieu> > ey; 
t.back() .push back(ey); 
for(z=0;z<nbZ;z++) 
{ 

ArRef<Milieu> m= = new_Milieu(); 

tback().back().pushback(m); 

m->setPosition(x ? y,z); 

> 

> 

> 
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for(x=0;x<nbX;x++) 
{ 

for(y=0;y<nbY;y++) 
{ 

for(z=0;z<nbZ;z++) 
{ 

ArRef<Interface> i; 
double diffusivity = 0.5; 
//double velocity - 1.0; 
if(x>0) 
{ 

i=new_Interface(); 
i->setSource(t[x] [y] [z]); 
i->setTarget(t[x-l][y][z]); 

i->setDiflftisivity(dififiisivity); 

i->setVelocity(0.0); 

} 

if(y>0) 
{ 

i=new_Interface(); 
i->setSource(t[x] [y] [z]); 
i->setTarget(t[x] [y~ 1 ] [z]); 

i->setDiffusivity(diffiisivity); 

i->setVelocity(0.0); 

} 

if(z>0) 
{ 

i=new_Interface(); 
i->setSource(t[x] [y] [z]); 
i->setTarget(t[x][y][z-l]); 

i->setDiffiisivity(diffusivity/l 000); 

i->setVelocity(0.0); 

} 

} 

} 

} 



S31 



cerr « ArClass::find("Interface")">getNbInstances(true) 
« " interfaces" « endl; 
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S32 


cerr « ArClassrifindC'Milieu'O^getNblnstancesCtrue) 
« " milieux" « endl; 




S41 


ArRef<Viewer3D> v=NEW_MyViewer(); 

ArRef<B oundingB ox3 D> bbox=new_BoundingBox3D(); 

v->getScene()->computeBoundingBox(bbox); 

v~>setFarDistance(5.0*bbox->getMaxSize()); 

v->setNearDistance(0.00 1 *v->getFarDistance()); 

unsigned int w,h; 

Viewer3D: :getScreenSize(w,h); 

V">setWindowGeometry(w-400 ? 0 ? 400 ? 400); 

v->setMapped(true); 

return(sched); 
} 




S50 


int 

main(int argc, 
char ** argv) 

{ 

ArSystem3D init(argc,argv); 
re1xim(Mt.simulate(&simulationInit)); 

> 





5 
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Annexe X2 
Representation simplifiee 



Classe Interface 


Constructeur d'objets 


Methodes specifiques 


setVelocityO 




getVelocityO 




setDiffusivity() 




getDiffusivityQ 




setSource() 




getSource() 




accessSourceO 




setTarget() 




getTargetQ 




accessTarget() 




getOtherEndO 




accessOtherEndQ 


Methode activate 


activate() 




Classe Milieu 


Constructeur d'objets 


Methodes specifiques 


setDensityO 




getDensityO 




setColor() 




getColor() 




addSourceQ 




removeSource() 




addlnterface() 




removeInterface() 




getNblnterfaceQ 




getlnterface() 




accesslnterface() 


Methode activate 


activate() 
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REVENDICATIONS 

1 . Dispositif (D) de simulation du monde reel, propre a etre implante dans un ordinateur 
5 (C) agence pour supporter en mode multi-taches une programmation par objets 

actives, representatifs de systemes a simuler, caracterise en ce qu'il comprend un 
logiciel de simulation par objets de 1'evolution conjointe de certains au moins desdits 
objets actives, comportant : 

- des objets d'etat contenant chacun au moins une donnee d'espace et/ou de temps 
1 0 et/ou au moins vine donnee de propriete, definissant un etat courant, 

- des objets d'interaction contenant chacun la designation d'au moins Tun desdits 
obj ets d'etat et d' au moins une fonction applicable a Pun au moins de ces objets d'etat, 
et definissant a chaque instant la topologie du systeme simule, 

- un gestionnaire de simulation capable de travailler par sequences sur une selection 
1 5 d' objets d' interaction, et a activer chaque objet d' interaction une unique fois lors de 

chaque sequence, selon un ordre variant de fa9on au moins partiellement aleatoire 
d'une sequence a l'autre, de maniere a appliquer chacune de ses fonctions a 1'etat 
courant de chaque objet d'etat qu'il designe pour faire evoluer son etat vers un nouvel 
etat courant. 

20 

2. Dispositif selon la revendication 1 , caracterise en ce que ledit logiciel de simulation 
comprend des objets d'interaction interne, aptes a contenir chacun la designation 
d'un seiil objet d'etat et d'au moins une fonction applicable a cet objet d'etat, et des 
obj ets d ' interaction mutuelle, aptes a contenir chacun la designation d ' au moins deux 

2 5 objets d'etat et d'au moins une fonction applicable a des donnees de propriete de ces 

objets d'etat designes. 

3. Dispositif selon Tune des revendications 1 et 2, caracterise en ce que ledit logiciel 
de simulation est agence pour modifier certaines au moins desdites fonctions en 

3 0 fonction d'au moins une donnee de propriete d'au moins un objet d'etat associe. 
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4. Dispositif selon Tune des revendications 1 a 3 ? caracterise en ce que ledit logiciel de 
simulation est agence pour selectionner certaines au moins desdites fonctions en 
fonction d'au moins une donnee de propriete d'au moins un objet d'etat associe. 

5 5. Dispositif selon Tune des revendications 1 a4, caracterise en ce que certains au 
moins des objets d'etat comprennent une donnee de propriete representant une 
grandeur intensive. 

6. Dispositif selon Tune des revendications 1 a 5, caracterise en ce que certains au 
10 moins des objets d' interaction ont une fonction faisant intervenir une grandeur 

extensive ou intensive. 

7. Dispositif selon Tune des revendications precedentes, caracterise en ce que certains 
au moins des objets d'etat comprennent des sous-objets d'etat. 

15 

8 . Dispositif selon la revendication 7, caracterise en ce que certains au moins des objets 
d'etat comprennent des sous-objets d'interaction operant sur lesdits sous-objets 
d'etat. 

2 0 9. Dispositif selon l'une des revendications precedentes, caracterise en ce que ledit 
logiciel de simulation comporte des classes d'objets definissant des structures 
d'objets d'etat et d'objets d'interaction, lesdits objets d'etats et d'interaction etant 
derives de ces classes par instanciation. 



25 10. Dispositif selon l'une des revendications precedentes, caracterise en ce que ledit 
logiciel de simulation comporte un ordonnanceur propre a fonctionner selon l'un de 
deux modes choisis parmi un mode en temps reel, dans lequel il fonctionne selon une 
frequence choisie, et un mode en temps virtuel, dans lequel il fonctionne de fafon 
periodique mais pendant des durees variables d'une periode a Pautre. 
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Debut sequence n 



Tirage au sort d'un objet 
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Activation et execution 
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NON 




Fin de la sequence 







n = n+ 1 
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